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ABSTRACT 


A high resolution, multi-level, primitive equation (PE) model of the California 
coastal region is initialized with a temperature field analyzed from real data collected 
during the OPTOMAII cruises of June and July,1984, for the purpose of forecasting the 
movement of thermal features in the region. The results are compared to the observa- 
tions and to the forecast experiments of Rienecker er al.(1987), since they initialized their 
quasi-geostrophic (QG) model with the same OPTOMAII1 temperature data. Key pre- 
diction features include an anticyclone and cyclone pair, and an offshore “jet” that was 
formed between the pair with velocities on the order of 60 cm/sec in the upper ocean 
region (< 225m). The temperature front associated with this “jet” 1s traced at the 85m 
level, in a time series from day 0 to day 14, as is the perturbation pressure field. 
Translational velocities of this frontal feature are on the order of 5-10 cm/sec in a 
southward direction, which 1s consistent with those observed. Some quantitative differ- 
ences between the PE model prediction and the QG model are found. Based on these 
results, it is feasible that frontal movement in the California coastal region can be fore- 
cast by a multi-level, high resolution PE model, given synoptic data for initialization. 
However, manv more studies are needed to understand the dynamics and robustness of 


the present model predictions. 
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I. INTRODUCTION 


A. DESCRIPTION OF THE CALIFORNIA CURRENT SYSTEM (CCS) 

The California Current is more accurately described as a system because of its 
broad, filamentous and variable nature. It is essentially the equatorward leg of the 
North Pacific gyre which 1s defined by the strong Kuroshio Current on the western 
boundary, yet the California Current 1s much less well-defined, as are most true eastern 
boundary currents. Broad and slow in the mean, its southward trend 1s often interrupted 
bv offshore “squirts” and meanders which form cold and warm eddies, often in dipole 
pairs. Chelton (1984) states that the core of the mean California Current 1s located 
around 100-200 km offshore in the upper 200 meters, with velocities ranging from 5 to 
14 cm/sec depending upon season and location. These velocities are computed 
geostrophically using 500m as a reference level, 1.e., the depth at which the current 1s 
assumed to be zero, and using over 30 years of hydrographic survey data from California 
Cooperative Oceanic Fisheries Investigations (CalCOFI) in the CCS between San 
Francisco and Point Conception. Velocities up to 40 cm/sec in the equatorward surface 
current are found further north in the Cape Mendocino region in the summertime, and 
the current can penetrate as deep as 300-400 meters. The prevailing southward winds 
force the surface waters offshore due to Ekman drift in the summer season which results 
in upwelling at places such as Point Sur and Point Arena, bringing up nutrient-rich, 
colder water from greater depths. 

There also has been noted (Hickey,1979) a poleward undercurrent called the 
California Undercurrent, with mean velocities of about 5 to 20 cm/sec in summer which 
hes below and inshore of the southward-flowing California Current. This current can 
surface during winter months when winds are not controlled by the dominant high 
pressure of summer over the North Pacific, and it is then known as the Davidson Cur- 
rent. 

The CCS is important to the climate of the western coastal states. Although not 
as Well-defined as western boundary currents such as the Gulf Stream, and certainly not 
as deep, it can bring up nutrients and cold water from the depths during an upwelling 
event, it can provide “air conditioning” to the beaches from onshore breezes, and it can 
provide moisture for the surface boundary layer, assisting cloud formation and fog. 


Fronts associated with temperature, salinity or density gradients form at the boundaries 


of the CCS, affecting acoustics performance and the fisheries. The prediction of frontal 
formation and movement in the California coastal region becomes therefore, an impor- 


tant goal of ocean modelers. 


B. HISTORY OF OCEAN PREDICTION IN THE CCS 

Most ocean modelers have used seasonal climatology as input data to their models, 
but the model simulations in these cases cannot be considered ocean predictions. A 
good example of such a study is that of Ikeda and Emery (1984), who showed that a 
particular baroclinic jet, typical of that which is observed in the CCS north of Cape 
Mendocino in the summertime, 1s baroclinically unstable. Rienecker ef al/.(1987) used 
observed data to initialize their quasi-geostrophic (QG) ocean model to create a dy- 
namical interpolation and forecast of mesoscale features in the CCS. These initial data 
Were taken from OPTOMAII, the Ocean Prediction Through Observation, Modeling, 
and Analysis program, which was designed for that very purpose. Six cruises in the same 
area approximately 180 kilometers offshore from Point Arena were conducted during 
June and July 1984, each collecting data with XBT drops and CTD casts. Figure | de- 
picts the cruise area of one of the legs of OPTOMAII. These cruises were timed to 
provide synoptic data for initialization and verification of model forecasts. Rienecker 
et al{1987) were able to make successful dynamic forecasts using a six-level QG model, 
but the success for some of the time periods required the inclusion of wind stress curl 


forcing in the model. 


C. THESIS OBJECTIVE 

With the goal of applying a primitive equation (PE) model to the problem of coastal 
ocean prediction. I have done essentially the same thing as Rienecker ef al.(1987), using 
the data from OPTOMAII to initialize a PE model of the CCS developed by 
Haney(1985) and modified by Batteen(1988). From physical considerations and simple 
scale analysis, it is clear that a multi-level PE model with active thermodynamics has 
greater potential for coastal ocean predictions than does a QG model, and this is the first 
thesis to specifically apply such a model to this problem. However, to keep this first 
study as simple as possible, wind forcing and variable bottom topography are both 


omitted. 
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Figure 1. Selected Stations from OPTOMIAI1, Leg DIT, (after Wittmann ez 
al..1985) 


li. THE MODEL 


A. DESCRIPTION 
The numerical model is a ten-level, primitive equation model of a baroclinic ocean 
on an f-plane, based on the hydrostatic, Boussinesq, and rigid lid approximations. The 


governing equations, Written in sigma coordinates and standard notation, are as follows: 
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For the finite differencing, a space-staggered B-scheme patterned after Arakawa and 
Lamb (1977) is used in the horizontal, and a sigma coordinate system defines the verti- 
cal. The boundary conditions at the sea surface consist of no heat flux and no wind 
stress. At the bottom of the ocean there is also no heat flux, and the stress is computed 
from the currents in the bottom level assuming a quadratic drag law with an inflow angle 
of 10 degrees and a drag coefficient of 1.2 x 10-3 (non-dimensional). 

The eastern boundary of the model domain represents a straight idealization of the 
continental shelf of the western United States with a zero-velocity condition imposed, 
while the northern, southern and western borders are open boundaries which use a 
modified version of the radiation boundary conditions of Camerlengo and O’Brien 


(1980). The total depth of the ocean model is a constant (3800m), with the midpoints 


of the ten levels located at depths ranging from 10 to 3100m (see Table 1). Figure 2 
depicts the actual topography in the model domain, although a flat bottom condition 1s 
used in the model. The use of variable topography, and an examination of its effects 


on coastal ocean predictions, 1s considered to be a separate study. 


B. TEMPERATURE DATA AND RESOLUTION 

Development of the temperature input field was done by subjective analysis and in- 
terpolation of data from OPTOMAII, Leg DII, which was conducted from 30 June to 
10 July, 1984. Seventy-four stations were chosen which were taken between 30 June to 
5 July, producing a quasi-synoptic pattern. The data domain (Figure 3) was a rectangle 
enclosed within latitudes 39.09 and 37.07N and longitudes 124.16 and 127.02W, which 
covers an area of 31,000 square kilometers, or 203 by 155 km. Final gridpoint resolution 
in the model was 3.23 km for delta-y and 4.23 km for delta-x, and the data field was 
framed by borders of 34 km wide on the zonal boundaries and 26 km wide on the 
meridional boundaries. Final size of the model domain (original data area plus borders) 
was 271 by 207 km, with the borders being used to blend in the temperature anomalies 
in various manners which will be discussed later in this section. Each level of the model 
was represented by 65 by 65 gridpoints. Since observed XBT and CTD data are limited 
to the upper 450 meters with only a few stations reaching deeper than that (to 735m 
depth), the bottom four model levels were constructed from the mean and the anomalies 
at the 450-meter level. 

After the observed temperature fields in the area of data, Figure 1, were hand- 
contoured for the upper six model levels (1.e., those at depths shallower than 450m), a 
13 by 13 gridpoint matrix was constructed from each of the resultant contour fields. Bv 
using linear interpolation, these matrices were refined to 49 by 49 gridpoint matrices, 
representing the inside area of Figure 3. The north, south, and western boundaries of 
the PE forecast model were treated as open boundaries in the manner described by 
Camerlengo and O’Brien (1980), as stated earlier. To reduce the gradients in the vicinity 
of these open boundaries (this is not a study of boundary forcing), the temperature 
anomalies (departure from the average at a given level) were extrapolated from their 
values at the outer edge of the data region (Figure 3) to one-half that value at the outer 
margin (model open boundaries). 

The eastern boundary of the model was treated differently since it represents a 
zero-velocity coastal margin. An isothermal condition at each level was desired on this 


boundary, in order to ensure that the across-shore thermal wind was zero on the coastal 





Figure 2. Bathvmetry of study domain: (The rectangle represents the PE model 


domain) 


margin. To accomplish this, the alongshore mean of the observed temperatures was 
calculated at the shoreward edge of the data region (Figure 1) and its across-shore slope 
was used to extrapolate the mean temperature to the eastern boundary which was lo- 
cated 34km from the edge of the data. The rest of the border was filled in by interpo- 
lation between the values at the edge of the data and this extrapolated temperature at 
the coastal boundary. This procedure was selected from among several possible methods 
of filling in the data between the coastal boundary and the shoreward edge of the 
OPTOMAT! observations. The result of applying this analvsis procedure at the $5m 


level 1s shown in Figure 3. The area of observed data is represented by the inside frame 


in this depiction, with the temperature contours superimposed. The contouring outside 
the framing represents the extrapolation of the data as described above. The main 
thermal feature at this level 1s the two-degree temperature gradient in the central- 
southern portion of the domain, separating the warmer northern portion from the cooler 
southern portion. This strong thermal gradient 1s accompanied by a rather intense 
geostrophic flow, or jet, directed offshore in the center of the domain. 

The bottom levels of the model had to be initialized from temperature extrapolations 
from the deepest model level that contained ample information to make an analysis. 
The mean temperature of this deepest level, 450 meters, was calculated and extrapolated 
downward exponentially with a decay scale of 700 meters, to a depth of 3100 meters. 
The temperature anomaly at this same depth of 450 meters was also extrapolated 
downward, but using a Gaussian formulation to virtually eliminate any measurable 
anomaly below 1200 meters. Table | displays the calculated temperature means at the 
various levels and compares these actual values to climatological values that have been 


used in the absence of real data. Also displayed are the temperatures in the warm eddv 


Table 1. COMPARISON OF TEMPERATURE (C) 
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noted at point A on Figure 3. The climatological temperatures were used in a modeling 
study bv Batteen (1988), and show a profile decreasing from 15 degrees at the surface to 


2 degrees at 4500 meters in an exponential fashion, with no anomaly at all. The actual 
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Figure 3. Model Domain Temperature at 85-Meter Depth: (inside frame repres- 


ents data area) Contour interval is 0.5C. 


field in this experiment ranged from 14.4 degrees at the 10-meter level to 2.09 degrees 
at 3100 meters. The major difference is of course, the anomolous temperatures extend- 
ing to over 1200 meters. The temperature gradients associated with such anomalies de- 


termine the initial velocity fields through the thermal wind relation. 


Figure 4 depicts the initial vertical temperature profile in the warm eddy (point A in 
Figure 3) and the mean temperature profile of the whole domain. The mean temperature 
profile of the entire domain decreases with depth, from 13.6 C at the 10m level to 2.09 
C at 3100m, while the profile taken at point A of Figure 3 is representative of the warm 
eddy. At 10m the temperature is 14.45 C, and it also decreases with depth to 2.09 C at 
3100m, but the warmer temperatures persist to 1200m or deeper. Figure 5 displays this 
more Clearly, with the upper 40m being a full degree or more warmer than the mean and 
the difference at 1200m only about 0.3 C. These figures show that the vertical extrapo- 
lation procedure was effective in spreading the temperature fields (means and anomalies) 
With depth. In the future, more sophisticated methods such as dynamical or empirical 
modes may be used for this purpose. These two figures clearly demonstrate the striking 
differences between using a climatological temperature profile as an initial condition, and 


an actual temperature field. 
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Figure 4. Temperature Profile: Warm eddy (Point A) vs Mean profile for entire 


domain 
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Il. MODEL RESULTS 


A. INITIAL TEMPERATURE, CURRENT AND PRESSURE FIELDS 

The quasi-svnoptic temperature field of 30 June-5 July has a central date of 3 July 
which 1s considered to be day 0, or the initial start time for this model simulation. Day 
6 then corresponds to 9 July and day 14, to 17 July. Figure 6 shows the horizontal 
temperature fields at day 0 at depths of 10, 85, 150 and 225 meters. While each level 
clearly shows its own distinct thermal signature, there 1s considerable vertical continuity 
as well. At each level a warm mass or wedge in the north-central portion of the domain 
1S In contact with a colder section in the southern portion of the domain, producing a 
gradient, or thermal front, between the two masses. The time evolution of this gradient 
will be described in detail in Section III.B since it is readily identifiable and since it 1s 
ultimately responsible for the strong initial currents in the upper ocean. These currents, 
computed from the thermal wind relation and the assumption that the verticallv aver- 
aged current 1s zero, are depicted in Figure 7 at the 10-meter level. The solid contours 
represent onshore zonal velocity and poleward alongshore velocity. The strong offshore 
jet in the center of the domain has typical speeds greater than 40 cm/sec with a maxi- 
mum value of just over 70 cm/sec. Rienecker er a/. (1987) also found offshore velocities 
of 60 cm’sec in the same region. The alongshore currents, southward in the eastern part 
of the domain and northward in the western part, are somewhat weaker with maximum 
magnitudes near 40 cm'sec. 

The actual flow is the resultant of the two velocity components and 1s more readily 
inferred from the bottom diagram in Figure 7, which shows the perturbation pressure 
field (dynamic height in centimeters) of the top level in the model relative to 2400m. In 
the diagnostic computation and display of dynamic height fields, 2400m represents the 
level of no motion, and will be discussed in Section III.C. The solid contouring in the 
northern part of the domain in this “P” diagram 1s indicative of a warm high pressure 
center with associated anticyclonic flow. The dashed contours in the lower half repre- 
sent lower pressure and thus, cyclonic flow around the low centers. This depiction 
clearly reveals what happens to produce a “squirt”, or offshore jet. The anticyclone and 
cyclone combine to act as a “nozzle”, accelerating the water seaward at high velocities. 
The “U” diagram shows the strongest offshore speeds in the center of the domain, which 


corresponds to the area of strongest gradient on the “P” diagram, between the two 
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Figure 6. Input temperature fields at 10, 85, 
interval(ci)=0.5 C. 


130 and 225 meters: contour 


pressure centers. In contrast, the maximum onshore velocity is onlv 40 cm’‘sec, and that 


is generated in the southern border region which was imposed by the assumed damping 
of the temperature anomalies toward the boundaries. 
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contours represent high pressurc) 


Figure 8 shows sections of the velocity fields at day 0. The top figure shows the 
shoreward (U) velocity field in a north-south section through the center of the high 
pressure area While the bottom represents the alongshore (V) velocity, in a zonal section 
approximately halfway through the domain. The solid contours represent onshore flow 
(LU >) in the top diagram and poleward flow (V>0) in the bottom diagram. Taken 
together the two sections show the large-scale current patterns in the upper 2km of the 


model. These include the strong offshore jet (L < Ojnear v= $84km with weaker onshore 


14 


flow (LU >Q) both to the north and south of the yet. It can be seen that the large-scale 
currents associated with the initial temperature fields are significant primarily in the up- 


per 100Um. 
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Figure 8. Sections of Shoreward(U) and Alongshore(V) Velocity, Day 0: Taken 
through near center of domain, (ci= 10cm‘sec). (Depth scale shortened 


for easier compurison) 
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The above “snapshot” of the fields at day 0 (3 July) depicts the velocities and flow 
patterns produced solely by the temperature input and the PE (geostrophic) computa- 
tion in the model, at initialization. What follows is an examination of the evolution of 
these features with time, as observed and as predicted by the PE model. 


B. TIME EVOLUTION OF THE FIELDS 

A good overview of the time development of the flow field predicted by the PE 
model for the two weeks following the initial time is seen in the near surface temperature 
field (Figure 9). The thermal signature present in the main upper ocean front is distinct 
enough to recognize from day to day, so its movement can be easily seen. The 85m level 
is the best level to use for this purpose because of the strength of the gradient, although 
a strong gradient is also present at deeper depths (450m). In the center of the domain, 
this thermal gradient moved southward about 25-30 km in the first four davs, indicating 
a translation speed of 6-7 cm/sec for the front. What is even more evident 1s the 
intensification of the gradient, 1.e., from an initial 2-degree change in 30 km, the gradient 
has tightened to 2 degrees in 15 km, or doubled in intensity. The portion of the front 
nearer the coast has pushed northward about 80 km in the first 6 days, while the rest 
of the front has moved south and west, stretching and thus, tightening the gradient. 
During this two-week period that the front intensified, small mesoscale features or Warm 
cells were seen on the warm side of the front traveling from the northeast part of the 
domain southwestward, then northwestward along the front. At the same time the cold 
water to the south of the front pushed eastward and then northeastward, reorienting the 
front in a more north-south direction. 

A comparison of the time evolution of the perturbation pressure at 85m (left column 
of frames) with that at 225m (right column) is shown in Figure 10 from day 0 to day 14. 
At 85m, day 0, the dynamic height is a maximum of 22.5 cm in the center of the high 
pressure area, and a minimum of -15 cm at the center of the low near the eastern 
boundary. The total difference is 37.5 cm, while the 225m-level range 1s from 17.5 to 
-12.5 cm, or a 30-cm difference. By day 14, the difference has increased slightly to 40 
cm at 85m and 32.5 cm at 225m, both increases due to an intensification of the low 
pressure center. 

In those two weeks, the center of the warm high pressure at the 85m level moved 
60 km to the south while the axis of the jet between the warm high and the cold low has 
“backed” from westward to southwestward, paralleling the coast. This is due to a 


northward migration of the low center relative to the high center, keeping the “jet” 


16 


(us) @2U0) 619 





(ws) 82U0]81Q 


(us4) 82UD] 81g 


(ws4) aU] 81Q 


121 $1 


131 


Distances (km) 


181 151 {21 91 


214 


271 261 


Distonce (km) 


5 C) 


V intervals, ci= 


days 0-14(two-day 


Sm Temperature front, 


8 


Figure 9, 


a 








between them. Consequently, the zonal component of the jet is weakening at day 14, 
while the meridional component is strengthening. The 225m level dynamic height field 
is less intense but the pattern is very similar to the upper levels. The southward 
displacement of the jet is about 5 km less, suggesting a north-south tilt with depth. This 
can be seen also in Figure 11, which shows the north-south section of the U-field and 
the east-west section of the V-field, both taken at day 5. Figure 12 shows the same ve- 
locity sections, only at day 14. It can be seen that at day 5, the westward zonal jet is 
located further south at the surface than at 800m depth. Thus, the jet is tilting towards 
the warm water with increasing depth. Such a tilt is suggestive of baroclinic instability, 
but a thorough analysis of the actual energy conversions is required to investigate this 
possibilitv. Regardless of the mechanism, the PE model predicts a strong development 
of the thermal field and associated currents during the two weeks following initialization. 
According to the model, the cold low pressure system located south of the main jet in- 
tensifies and moves east-northeastward, reorienting the jet to a more north-south direc- 


tion. 


C. VERIFICATION/COMPARISON WITH QG MODEL 
1. Domain and Initialization 

Before making comparisons with Rienecker ert al. (1987), the domains of each 
model run must be matched geographically in order to have meaningful comparisons. 
Figure 13 shows the data domain used in this paper (solid lines) and the domain used 
by Rienecker er al. (1987) (dashed lines). (See also Figures 1,2 and 3 of this paper.) 
Those authors used data from OPTOMAII to initialize a 150km square domain for their 
QG model while the PE model used a rectangular domain, including all of the data out 
to longitude 127.02W. These dimensions do not include the border regions of either 
model. 

Figure 14 shows this relationship of orientation on the sea _ surface 
temperature(SST) and surface dynamic height(SDH) at initialization. Figures 14a and 
14c are adapted from Figure 3 of Rienecker et a/.(1987). Figure 14a depicts a sea surface 
temperature map made by objective analysis using data from two cruises of 
OPTOMAI]I, covering an asynoptic time frame of 23 June to 10 July, 1984. The framed 
area represents the data domain of the PE model, which used data taken 30 June to 5 
July, 1984, OPTOMAII, Leg DIT only, for a quasi-synoptic picture. Figure 14b shows 


that same area, with the temperature contours done initially by subjective analysis and 
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Figure 11. Sections of Shoreward(U) and Alongshore(V) Velocity, Day 5: Taken 


through near center of domain, (c1= 1Ocm/sec). (Depth scale shortened 
for easier Companison) 


plotted for the 10m level. The similarity in the patterns is evident even though separated 


by 10 ineters and prepared by different techniques. 
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Figure 12. Sections of Shoreward and Alongshore Velocity, Day 14: Taken 
through near center of domain, (ci= !0cm,sec). (Depth scale shortened 


for easier comparison) 


Figure I4c depicts the surface dvnamic height(SDH) relative to 450m. with the 


asvnoptic tume frame as in Figure !4a, as computed from observations by objective 


=e See a Dy VE is 
= i \\ f 
= Seeley =z Amn coon Mendocino 
WS TICS 





Fieure 13. Nlodel Domains( PE = solid, QG = dashed) 


analvsis. Figure I4d shows the perturbation pressure at 10 meters relative to 450m 
computed from the PE model. Again. the similar patterns are evident despite the slight 
differences in svnopticity of input data. These figures provide evidence that both the 
QG and the PE models were initialized with essentially the same data; therefore the re- 
sults of each should be a good comparison of the dynamics involved, despite the slight 


differences in synopucity. 
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2. Verification 

Rienecker er a/l.(1987) constructed SDH for each of the six legs or cruises from 
OPTOMAII (total time span 5 June to 5 August, 1984) to compare with the QG model 
runs. Three of these will be borrowed from their Figure 5 for use in this paper. Figure 
I5a depicts the observed SDH from Leg AIII, 5-13 July, which corresponds to day 6 of 
the PE model run. This field was computed by objective analysis, Rienecker er al. 
(1987). Figure 15b shows the PE model-predicted surface P-field relative to 450m of day 
6 aS a comparison. Figure I5c shows the 50m (non-dimensional) dynamic 
streamfunction, relative to 750m, as predicted for day 6 (their day 193) by the QG model 
of Rienecker er a/. (1987) using wind stress curl forcing estimated from FNOC wind an- 
alyses. Figure 15d shows the P-field at 50m, also relative to 750m, predicted for day 6 
bv the PE model of this paper. The inset squares on the two right-hand figures represent 
the area of the PE model domain that corresponds to the 150km sq area of the observed 
and QG model prediction of the left figures. On both PE figures, the center of the high 
pressure is in the upper left corner of the inset box, which matches fairly well with the 
actual location of the high on Figure 15a, but the center of the low pressure area is 
about 30km east of the actual position on the verification panel. This difference in 
position of the cyclonic center effectively removes the “crimp” out of the “jet” that 1s so 
apparent in Figure I5a, and consequently the PE pressure gradient runs diagonally 
rather than zonally. Figure I5c, the QG-predicted streamfunction at 50m, shows more 
of the zonal flow at day 6 than does the PE P-field. The QG-predicted “jet” orientation 
here is similar to the verification panel, significantly better than the PE pressure field of 
Figure 15d, but the location and axis of the cyclone-anticyclone pair predicted by the 
PE pressure field is nearer reality than the QG streamfunction predicts. The SDH ver- 
ification field, Figure 15a, shows the “jet” between the high and low centers to have 
moved further offshore than the PE model predicted, and the feature itself distorting into 
a very tight gradient with an abrupt 90 degree direction change to the south, then an- 
other shift back to the west. This offshore motion 1s not predicted by the PE model be- 
cause the low pressure center has moved northward further in the model, causing a 
north-south “jet” axis shift between the high and low centers which eliminates much of 
the Zonal flow. 

Figure 16 is the same as Figure 15, except that prediction day 14 rather than 6 
is shown. The verification SDH in panel ”a” 1s taken from Figure 5 of Rienecker er al. 
(1987), and is based on the data collected on Leg P of OPTOMAII. This is the most 
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synoptic data of the program, having been taken by aircraft (P-3) on 18 July. 1984, in 
about 8 hours. Due to depth lnutations of data, a 30Um reference level was used to 
compute the SDH of the verification panel, so the same reference level was used for 
computing the day 14 surface P-field predicted by the PE model (panel b) for better 
comparison.(Note that dav 14 of the PE prediction is actually 17 July. one dav before 


the verification panel, but the elapsed time of model run is the same for PE and QG, 


tJ 
ai 


with one dav difference for the start time, so it 1s used for comparison.) Panels c 
(QG-predicted streamfunction at 50m) and d (PE predicted P’ at 50m) both use a 750m 
reference level, however. 

The high pressure center is still located in the northwest corner of the inset of 
the PE prediction panel, but the center of the high in the verification has moved west, 
out of the panel domain. The low center or cyclone is predicted by the PE model to be 
about 30km further south than the verification shows. (Recall from Figure 15a and b that 
the dav 6 prediction was north of the verification). The gradient or “jet” is nearly 
bisecting the domain meridionally on the prediction, while the verification still shows 
zonal flow in the northern third of the domain. The intensity of the predicted cyclone 
in Figure 16b has pulled the gradient tightly in and away from the anticyclone while the 
verification in Figure 16a shows the gradient to have weakened since day 6. The im- 
provement in synopticity of the later data (8 hours vice 8 days) may have something to 
do with this apparent weakening of the analyzed gradient. 

The QG prediction at 50m, shown in Figure 16c, again was done with wind 
stress curl forcing added. The gradient runs diagonally through the domain in the gen- 
eral direction of the verification field, northeast to southwest, but without the cyclonic 
and anticyclonic curvature. The 50m PE prediction is much the same as the surface 
prediction, with the gradient running basically north-south, but it finishes in the south 
with the anticyclonic turn displaved by the verification SDH of panel a. It appears that 
a combination of QG and PE predictions would come closest to the verification on day 
14. 

3. Discussion and Interpretations 

Rienecker er a/.(1987) chose a reference depth of 750m for a level of no motion, 
upon which to base the dynamic heights calculated and predicted by their QG model. 
In contrast, the PE model does not assume a particular reference level, but rather as- 
sumes that the vertically averaged current is zero. According to the model currents in 
Figures 11 and 12, this results in an approximate level of no motion near 1200 meters. 
Figure 17 shows the perturbation pressure field at 450 meters at initialization (day 0) of 
the PE model, and the resultant U velocity field. The vertical continuity is apparent with 
the persistence of the cyclonic-anticyclonic features which are seen at shallower levels, 
still present at 450 meters, albeit in weaker form. The downward extrapolation of the 
450-meter temperature anomaly is responsible for this dynamic picture present at this 


depth, and the dynamic height and, consequently, velocities, could be altered by using 
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different vertical extrapolation techniques. Figures 4 and 5 show a warm-anomaly of 
nearly one degree at 450m, however, and to terminate this anomaly abruptly, which 1s 
really closer to two degrees total differential across the front. does not seem to be a 
physically correct procedure. Without observations below 450m, it is not possible to 
know what the real temperature anomalies at great depths are in any given Case. 


Whether anv of the differences between the QG and the PE model results can be 


attributed to this difference in level of no motion cannot be determined without further 
sensitivity studies, as will be discussed in Section IV.B. 

For the initial conditions used in this study, the PE model was sensitive to 
changes in the open boundary conditions, and several variations were tried unsuccess- 
fully before a satisfactory condition was found. Virtually all the conditions chosen ex- 
cept the one that worked became unstable very quickly, in two or three days model run 
time, and produced poleward coastal jets of 500 cm/sec or more. Appendix B lists the 
open boundary conditions tried in detail, since this paper is not a study of open bound- 
arv conditions. Without further sensitivity tests it is not possible to know whether one 
of the open boundarv condition formulations is unstable (unlikely) or whether the linear 
computational stability condition Cores <1) was violated by the jet, or whether some 
other instability occurred. . 

The isothermal condition imposed at each model level along the coastline re- 
sulted in varving degrees of slope of the isotherms between the boundary of the data field 
and the coastline. These slopes produced thermal wind inflow from both the north and 
from the south (e.g., Figure 3). Since the slope of the isotherm in the south was greater 
than in the north, inflow from the southern boundary tended to dominate in this border 
region, simulating an intense poleward coastal jet that was subsequently predicted to 
intensifv with time. A different procedure at the eastern boundary would have produced 
a varving amount of flow from the south or from the north in the border region, result- 
ing in a different model initial condition and perhaps a different model prediction. An- 
other method of treating the eastern boundary 1s to “open” it, like the other boundaries, 
with no coastline. The robustness of the model results to the open boundary conditions 
and to the method of extrapolating the data to the various model boundaries needs to 


to examined further, as discussed in Section IV.B. 
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IV. SUMMARY AND RECOMMENDATIONS 


A. SUMMARY 

This PE model run utilized observed temperature data gathered by Leg DII of 
OPTOMALII, 30 June-10 July, 1984 as initialization for the purpose of predicting ther- 
mal feature movement in the California coastal region. The specific area of coverage, a 
31,000 square kilometer domain 180 km off of Point Arena, was sampled in intervals 
condusive to verification of model predictions, and Leg AIII, 5-13 July, and Leg P of 
18 July were used as verification of the PE model predictions of days 6 and 14. Rienecker 
et al, (1987) also used OPTOMAI]I in the same manner for the initialization of their QG 
model, so their results were used extensively for comparison as Well. 

The observed initial offshore jet velocities of 60 cm/sec or more in the PE model 
echo the findings of Rienecker et al. (1987), whose model produced similar surface jets. 
These jets were verified by satellite images that clearly showed at least the presence of 
jets and meanders, if not the velocities. The input surface temperature fields and surface 
dvnamic height and pressure surfaces of each model were matched with objective 
analvsis-calulated verification panels and found to be basically equivalent at the outset 
of model simulations. The differences between the two models at initialization were 
reference level differences and model domain differences, all of which were accounted for 
in the comparisons. The QG model's dynamic streamfunction at 50m was compared 
with the pressure perturbation of the PE model at 50m. Both models were close to ver- 
ification, even though the QG model was driven by observed winds and the PE model 
was unforced. For example, both models were within a few kilometers on the 6-day 
prediction of the anticyclone, but the QG prediction of the cyclone was 40-50km south- 
west While the PE model put the cyclone 30km northeast of the verified position. Con- 
sequently, the jet axis formed by the gradient between the cyclone/anticyclone pair was 
shifted accordingly. Therefore, by day 6, axis angles between the cvclonic and 
anticyclonic centers were different between the two models which resulted from 
translational differences in the dynamic features, with the QG model predicting offshore 
movement while the PE model predicted southward and northward motion of frontal 
features. 

By day 14, the QG model 50m streamfunction represented the jet basically 


diagonally across the domain while the PE model showed a nearly north-south 
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orientation of the jet. The verification SDH of Leg P of 18 July showed mostly offshore 
flow in the northeast quadrant of the domain, shifting cyclonically to southward flow in 
the central portion, and turning anticyclonically back to offshore in the southwest 
quadrant. Based on these results and those from the time series of the 85m temperature 
field, it is feasible that frontal movement in the California coastal region can be antic- 
ipated and, with more sensitivity and case studies (below), could be predicted by a 


multi-level, high resolution PE model, given synoptic data for initialization. 


B. RECOMMENDATIONS 

This study represents but one simulation (hindcast) of a single synoptic case in the 
coastal ocean using a PE model. An important question that needs to be investigated 
further is how robust the simulated eddy evolution is to various modeling and data 
processing assumptions. For example, one would like to know how sensitive the simu- 
lated results are to the procedures used to extrapolate temperature data into the regions 
(borders, deep ocean) of no data. The eastern boundary is especially important in this 
particular synoptic case because there is a strong jet entering the domain from the 
shoreward edge of the data region. Where did this jet actually come from? Did it come 
down the coast from the north and turn anticyclonically into the data area or did 1t come 
up from the south and turn cyclonically into the region? The method of extrapolating 
the data toward the lateral boundaries, with particular attention to the eastern bound- 
ary, needs to be examined carefully and the sensitivity of the simulation to these proce- 
dures needs to be thoroughly investigated. The sensitivity of the model predictions to 
the method of vertically extrapolating the fields to the deep ocean should also be studied. 
Alternative methods of vertical extrapolation, such as fitting the anomalies to dynamical 
or empirical modes, should certainlv be investigated. 

Another important question to address is what are the dynamical mechanisms that 
are responsible for the evolution of the eddy fields in this synoptic case. In order to in- 
vestigate these mechanisms, a number of sensitivity studies and analyses should be car- 
ried out. An analysis of the energy and vorticity budgets would be useful to determine 
the relative importance of purely kinematic effects (horizontal redistribution of absolute 
vorticity) versus baroclinic development. The sensitivity of these various dynamical ef- 
fects to the particular data processing and modeling procedures should also be studied. 
Finally the role of wind forcing (important in the QG simulations) and variable bottom 


topography should also be systematically investigated. 
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It is therefore recommended that further PE model studies be conducted using ob- 
served data for initialization and verification of dynamic numerical predictions. The 
OPTOMA program is ideal for that purpose since data which were taken in a quasi- 
synoptic fashion are provided in intervals of 1-2 weeks, or, as in the case of Leg P of 
OPTOMALII, within eight hours. A PE model could take such data coverage for an 
operational area and predict the thermal conditions 1-2 weeks later, which would be a 


valuable planning tool for the conduct of Naval exercises. 
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APPENDIX A. OPTOMAI1 STATIONS, LEG DII, 30 JUNE-10 JULY 1984 


XBT HEADERS: 

STA # DATE TIME LAT LONG NPTS 
1 84183 JES: 36.47 -122.10 684 
Zs 84183 200 BOeo -l22.21 664 
&. 84183 300 SO059 ~iZ2s1 441 
4 84183 407 o/s -122.41 662 
5 84183 503 ei aus -Ie2?.50 682 
6 84183 603 o/ mS -122.58 6864 
i 84183 700 B/ al -P25566 684 
8 84183 ous S7920 -13 7 684 
3) 84183 906 Sess -lNeone> OOF 

10 84183 1014 Bee? len 3/ “O52 

Ws 84183 110s 37.42 -123.45 684 

12 84183 1243 of o8 123459 6864 

Ls 84183 1350 Ofea) ~124°10 684 

15 84183 AAs 7 38.09 -124.24 684 

16 84183 2036 woe -ies. 38 669 

17 84183 220 38.34 -124.45 668 

18 84184 19 38.45 -124.55 684 

19 84184 430 Jo.0¢ -12>.01 654 

20 84184 806 39.02 -125.04 684 

Zl 84184 1106 Sooo -I>.11 516 

ZZ 84184 1705 e3mo2, -125.19 684 

Zo 84184 1800 Homo -l25.21 684 

24 84184 1850 960.43 -125.23 684 

25 84184 ness)! Sono2 -125,7/ 684 

26 84184 2022 Soo -l2Z>.25 8671 

2/7 84184 2102 36220 -125528 684 

28 84184 2139 como -l75n28 §~654 

30 84184 2344 wowed =125.36 684 

oil 84185 228 S7250°=1252543 684 

i ee 84185 700 Si250 =l25a,21 667 

34 84185 800 oe ~l25. 13 OS4 

35 84185 900 37.40 -125.04 684 

36 84185 1000 37.46 -124.52 684 

a7 84185 diz 37.51 -124.40 684 

38 84185 1206 37.57 -124.28 684 

40 84185 1444 38.03 -124.30 684 

41 84185 1544 38.03 -124.44 680 

42 84185 1636 S605 -124.°59 684 

43 84185 1722 36205 -125.08 664 

44 84185 1810 BOe07 -125.21 684 

46 84185 2155 38.14 -125.54 684 

47 84185 2302 35,515 -126506 684 

48 84186 6 souce ~126.23 1679 

49 84186 102 No.2] “126.28 654 

51 84186 248 Zoe) =-126.32 664 

a2 84186 426 Sonos -126.15 684 

53 84186 514 38.43 -126.05 684 
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APPENDIX B. OPEN BOUNDARY CONDITIONS 
A. INTRODUCTION 


The eastern boundary of the PE model was closed because it represents the coastal 
margin. The open boundary conditions in question relate to the northern, southern and 
western borders of the model domain, but in particular to the southern boundary since 
strong inflow occurred there, causing rapid progression to instability. The changes to 
correct the instability problem were made consistently each time, to each boundary 
(other than the eastern one), but in this description only the southern boundary condi- 
tions on the northward velocity (V) are described to avoid repetition. To eliminate ter- 
minology problems, the following are stated for clarification: 

e Superscript “t” refers to time differencing at gridpoints 


e subscript “B” refers to space differencing between gridpoints, where “B” represents 
the open boundarv 


e “V" is meridional velocity component 
e “Cis northward phase speed computed near the southern boundary 
e VD 1s the velocity difference between adjacent gridpoints 


e CV is the computed phase speed used to determine inflow or outflow at a bound- 
ary, depending on sign(<0 = outflow, >0 = inflow) 


B. PROBLEM #1 

The initial model run developed a strong surface inflow jet from the south imme- 
diately, and went unstable (V> 100cm/sec) by day 4. Even to depths of 735 meters and 
below, velocities above 300 cm/sec were generated by day 10 of model run time. The 
boundary conditions on the meridional velocity component at the southern open 


boundary were determined by first computing the spatial difference, 
VD = Vig} — Vijh. (7 
In the event that VD = 0, the boundary was set according to 
Vg = Vea [8] 


If VD #0, the boundary was set by computing 
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Ay Cn Vi) 


ole i i 
so thatif CV < 0 (outflow), 
Vip = Vins [10] 
or if CV > 0 (inflow), 
eV een 


It was suspected that the instability might be related to the condition used when VD 
= 0 (Eqn.8), since it calls for an instantaneous equalizing of properties at the next grid 
position or location, or an automatic “inflow” condition. The first solution tried was to 
change the condition to essentially an “outflow” situation, when VD = 0. Thus 


Equation [8] was replaced by 


Vit) = Vous, [12] 


and the model was restarted at day 0 with the same temperature input. This change was 
designed to eliminate any excess “inflow” that might have caused the intense poleward 
jet and subsequent instability, but after four days of run time, the instability was still 
there and there was no noticeable change in results from the initial run, out to ten days. 

The second change then tried was in the same pattern, but changing back to an 


“inflow” condition when VD = 0. Thus, Equation [8] was changed to 
Ve = Vp, [13] 


and the model was again restarted at day 0. Again, after four days of run time, the in- 
Stability was still there, and there was no noticeable change in results from the initial run, 
out to ten days. It was therefore apparent that the instability observed was not due to 
the manner of handling “inflow” and “outflow” at the boundaries when VD = 0. Nev- 
ertheless, a third change was made in the “inflow’/"outflow” condition to improve the 
handling of the radiation borders. A combination of the first two changes was adopted, 


in the average form, 


Ve' = 0.5(Vga, + Vp) [14] 
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and while not solving the problem, it was felt to be a better form than the original [8]. 


C. PROBLEM #2 
The other possibility for the observed instability lav in the computation of “C”, the 
phase speed. The scheme used to this point for computing CV (i.e., eqn [7] and [9]) is 
rewritten here, 
t t—1 t-1 tI 
V4 a V4 Gi: (Va4y = V4) [15] 
At ~ Ay 
to show that it really looks like a forward time differencing scheme. The first solution 
tried was a backward time scheme between the time levels t and t-1. Thus, Equation [15] 


was replaced by 


Vee oat ae (Vig ee re) 
= —C lV ————_. 
At c Ay L16] 


Once again the model was restarted, with the latest change incorporated in all open 
boundaries, and once again, the velocity jet inflowing from the south became unstable 
in four days, with no noticeable improvement. 

Another solution was tried, which used a leapfrog scheme based at time zero of the 


following form: 


pith ptt tal 


[+] 
Baan” Bas Gi (Vien, + Vea vi | 





[17] 


2Ar Ay 2 | B+2 
This change worsened the instability that was observed in the first runs, and was rejected 
along with the other unsuccessful conditions, Equations [8], [12], [13], [15] and [16]. 

A final solution was tried, again using a backward time step between the time levels 


t+ 1 andt. Equation [15] was replaced by 


t+] +1 +1 
Vee = Vina —~_cy (Vea — V4) (18) 
Ar = Ay 


After this change of the boundary conditions, the model was restarted. This time the 
instability disappeared. The inflow jet from the south was still present but the magni- 
tude did not intensify as in the previous runs. Meridional velocity (V) peaked at 


40-50cm/sec while zonal velocity (U) reached 70 cm/sec at a grid location representing 


38 


the tightest gradient. This change became the run used for model comparison. More 
study needs to be done with these boundary conditions to determine the reasons for this 
instabilitv, and the recommendations section of this thesis suggests several kinds of ex- 


periments. 
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